Microscopic Calculation of the Dielectric Susceptibility 

Tensor for Coulomb Fluids 

L. Samaj 

Institute of Physics, Slovak Academy of Sciences, Diibravska cesta 9, 

842 28 Bratislava, Slovakia; e-mail: fyzimaes@savba.sk 

February 1, 2008 

Abstract 

In a Coulomb fluid confined to a domain V, the dielectric susceptibility tensor xv 
depends on the shape of V, even in the thermodynamic y — > oo limit. This paper 
deals with the classical two-dimensional one-component plasma formulated in an 
elliptic y-domain, equilibrium statistical mechanics is used. For the dimensionless 
coupling constant T = even positive integer, the mapping of the plasma onto a 
discrete one-dimensional anticommuting-field theory provides a new sum rule. This 
sum rule confirms the limiting value of xv predicted by macroscopic electrostatics 
and gives the finite-size correction term to xv- 
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1 Introduction 

Classical Coulomb systems are the prototypes for studying the effect of long-range inter- 
actions in equilibrium statistical mechanics. In dimension the Coulomb potential (ff at 
spatial position r— (r^, r^, . . . , r^), induced by a unit charge at the origin, is the solution 
of Poisson equation 

M%r) = -ej{r) (1) 
with £2 = 27r, 63 = 47r, etc. In particular, vav — 2 dimensions, one has 

<P%r) = -\n{\r\/ro) (2) 

where ro is a length scale, for simplicity, set to unity. A general Coulomb system con- 
sists of s pointlike species a — 1, . . . , s of the corresponding charges qa, embedded in a 
fixed uniform neutralizing background of density no and charge density po- The most 
studied one-component jellium or plasma (OCP) and the two-component plasma (TCP) 
correspond to s — 1 {qi — q), po — —qriQ ^ and to s — 2 {qi — —q2),Po — no — 0, 
respectively. The Coulomb system is confined to a domain V, which can be: 

(1) infinite, V R"; 

(2) finite, bounded by an impermeable hard wall (for the sake of simplicity, uncharged 
and with no image forces); 

(2) semi-infinite, i.e., bounded by a wall, but infinite in at least one of the parallel 
directions. 

The symbol {■■■)¥ will denote the canonical averaging over the domain V at the inverse 
temperature /3 = l/ksT, under the system neutrality condition. The microscopic total 
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number and charge densities at r are given by 

i 

M = '^i) (3b) 

i 

respectively, where the sums run over particle indices. The canonical average number 
and charge densities read 

nv{r) = {n{r))v, pv{r) = {p{r))v (4) 
At the two-particle level, one considers the two-body distribution 

nv{r, r^) = (^5(r- rj)6{r' - rk))v (5a) 
as well as its truncated form 



n 



'y(r, r') = nv{r, r^) - nv{r)nv{r') (5b) 



and the truncated charge-charge correlation function 

Sy{r\r^) = {p{r)pW))v " {p{r))v{pW))v (6) 
In the case of the OCP, S is expressible as follows 

Sv{r\r') = [n^(r, r') + ny(r)5(r- r')] (7) 

The long-range tail of the Coulomb force causes screening, and thus gives rise to exact 
constraints, sum rules, for the structure function S (see review |l[]). Their explicit form 
depends on the geometry of domain V. 

In the bulk regime, limy^pf," Sv{f\r') = S{\r— r'l) is known to obey the Stillinger- 
Lovett rules which imply the zeroth-moment (electroneutrality) condition 

J dr S{r) = (8) 



and the second-moment condition 

/3 



f3jdr {ryS{r) = J dr \r\^S{r) 



^ i=l,...,u (9) 



tt{u - 1) 

For the OCP, the fourth moment of 5* is related to the isothermal compressibility 0|,[01'IIH' 
so that the knowledge of the exact equation of state in dimension two provides the 
explicit form of the former [§,10. Very recently [|1^], the sixth moment of 5* for the 



2d OCP was derived using a renormalized Mayer expansion in density |TT|. The formal 
analogues of the fourth and sixth moments of S in the 2d OCP are, respectively, the 
zeroth and second moments of the truncated total number density correlation function 
{n{r)n{r')) — {n{r)) {n{r^)) in the 2d TCP, as was derived in Ref. from analogies with 



critical systems and in Ref. ||T3[ directly by using diagrammatic methods. 



For finite systems, the sum rule 

/ drSv{r\r^) = f dr' Sv{r\r') = (10) 
Jv Jv 

only tells us that the total charge in domain V is fixed. The information analogous to the 
second-moment formula (^ is contained in the dielectric susceptibility tensor xy, defined 

by 

x'^ = ^^{{p'nv-{nv{nv) (na) 

where 

= f drr'p(r) i = l,...,u (lib) 
Jv 

is the ith component of the total polarization in the system and \V\ the volume. Within 
linear-response theory, xv relates the average polarization to a constant applied field E, 
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(P*) = Yl'j=iXv^'' ■ With regard to ([10|), the tensor xv is expressible in two equivalent 

ways, 



Xv = J^(in J^dr2 r\riSv{ri\r2) 

l^dn l^dr2{rl-riySv{n\r2) (12) 



2\V\ 

As V ^ TV one would intuitively expect that, according to the bulk second- moment 
formula (|^), the diagonal components Xv = Xv = • • • 5 tend to the value 

xV^-f/d.M=sw = ^^ (13) 

However, this is not the case: due to surface effects, the xv limit depends on the shape 
of V. Its value is predicted by macroscopic electrostatics for homogeneously polarizable 
systems In the case of elhptic (z/ = 2) and ellipsoidal = 3) V^-domains, one 

introduces the depolarization tensor Ty 

V = 7 ^^^^ I dr-' r') (14) 



27r(z/ - 1) dr'-dri Jv 

where r is an arbitrary point in V. It is the fundamental property of elliptic and ellipsoidal 
domains that the tensor Ty is independent of point r eV, and depends only on the shape 
of V. With regard to Poisson equation (1), its diagonal elements Ty = Ty are constrained 
by J2i=i Ty = ei//[27r(z/ — 1)]. In the limit V IV, electrostatics yields 

= 27r(z/-l)Tf ^^^^ 
In the special case of a 2d disk or a 3d sphere, Ty is isotropic, so that Ty = e,^/[2nh'{h' — l)]. 
Consequently, 



in contradiction with the previously suggested estimate ([T3|). Formula (|T6D was checked 



and finite-size corrections were calculated in |14],[15| for the 2d OCP formulated on disk 



for the exactly solvable case of the dimensionless coupling T = Pq = 2 [jl6[, |T7 . 

This paper deals with the 2d OCP formulated in an elliptic domain, which includes 
the circularly symmetric disk and the limiting case, the strip. The statistics now depends 
on the only parameter - the coupling constant F (the particle density only scales appro- 
priately the distance). At F = even integer, the 2d OCP is mappable onto a discrete Id 
anticommuting-field theory |18|,[0. It is shown that sum rules come from specific unitary 



transformations of anticommuting variables, keeping the general form of the fermionic ac- 
tion. A nontrivial transformation of anticommuting variables is revealed to generate a 
new sum rule. For the elliptic domain, this sum rule confirms microscopically the asymp- 
totic formula (|1^) and gives the finite-size correction term to Xv explicitly in terms of 
boundary contributions. 

The paper is organized as follows. Section 2 recapitulates briefiy the mapping of 
the 2d OCP onto the Id fermionic model. Section 3 establishes the formalism of unitary 



transformations of anticommuting variables, which imply the known sum rules ||2^ . Com- 
plementary (to the author's knowledge as-yet-unknown) sum rules, nontrivial when some 
asymmetry of the V-domain is present, are given as well. In the key Section 4, using a 
special "nearest-neighbor" transformation of anticommuting variables, one derives a new 
sum rule providing a proper split of Xv ^'^^'^ asymptotic ([T5|) and finite-size correction 
parts. In Appendix, by explicit calculations in the 2d OCP on disk at F = 2, a test of the 
results is presented. 
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2 Mapping onto Id fermions 

The model under consideration is the 2d OCP of particles confined to a domain V. For 
point rEV, the cartesian {x,y), complex {z,z) or polar (r, 0) coordinate representations 
will be suitably used. The neutralizing background of density Uq = N/\V\ induces the 
one-particle potential 

0,(r)= / dV0^(|r-r'|) (17) 

For the elliptic \^-domain of interest, in the reference frame defined by the axis of the 
ellipse, x^/a? + y^/h'^ = 1, both tensors xv and Ty are diagonal. The fundamental 
independence of the depolarization tensor Ty of point r E V and the invariance of 
r) with respect to the refiection along the x oi y axis imply 

</)6(r) = const - ttT^^x^ - nT^y'^ (18a) 

with Ty = b/ {a + b),Ty = a/ (a + 6). The corresponding electric field reads 

E,{r) = - V0fc(r) = 27rT^xx + 27rT^^2/y (18b) 

where x, y are the perpendicular unit vectors in x, y directions. In the circularly symmetric 
case of disk a = b = R (radius), one has Ty = Ty = 1/2, so that 

<j)b{r) = -7rrV2, Eh{r) = vrr (19) 

The total Boltzmann factor associated with particle configuration {rj} is written as 



exp 



rnoJ2Mn)-Tj2(f)%\ri-r, 



i<j 



(20) 



For r = 2j, 7 being positive integer, it was shown in [|T^ that the canonical partition 
function of the system, Zy (we will omit in the notation the dependence on A^), can be 



expressed in terms of Grassmann variables = l)---)7)) defined on sites 

2 = 0, 1,...,A^— lofa discrete chain and satisfying the ordinary anticommuting algebra 
and integration rules as follows: 

Zv = Jv^pV^ exp[Svi^,^)] (21a) 

7(Af-l) 

Sv{^,i^) = E (21b) 
i,j=0 

Here, VipV^ = Hilo^ '^'^t^^ ■ ■ ■ d^'i^^d^i''^^ • • • d^i^'' aiid the fermionic action Sy involves pair 
interactions of "composite" variables 

N-l 

= E tf^.-el:^ (22a) 



+ + i) 
7V-1 



VI/, = i^^.-.i^^ (22b) 



Jl---i7=0 



i.e., the products of 7 anticommuting- field components with a given sum of site indices. 
The interaction strength is given by 



Wij 



I d^z z'z^w{z, z) (23) 

where w{r) = exp[rr;,o<^fe(r")]. Denoting by 

(...) = ^ I Vi,Vi . . . exp[Svi^, ^)] (24) 

the averaging over the Id fermionic system, the particle-number density (^) can be ex- 
pressed as 

7(7V-1) 

nv{r)=w{z,z) E (^^*.) (25) 

i,j=0 



the two-body distribution ( |5a| ) and its truncation ( |5bD as 



7(7V-1) 

riy(n,'^2) = w{zi,Zi)w{z2,Z2) E (^n%^i2^i2) -^I'^f 4'^2' (26a) 

niij2i2=o 



7(JV-1) 

nv{ri,r2) = w{zi, Zi)w{z2, Z2) ^ {E^^^ j^Ei^^ ,y z\' zf' z^^i' (26b) 

hjii2j2=0 

where (Si.^Hi^^'j^)^ = (Si^^Si^^^^) - For disk (|l|), since the 

Boltzmann weight w{r) possesses the circular symmetry, the interaction matrix Wij takes 
the diagonal form, 

Wij = 6ijWi, Jy '^'^^ r'^'uj{r) (27) 

Owing to the "diagonalization" of the action, Sy = J2i '^iWi'^i, only fermionic correlators 
(Hj^^jjSj2^j2 . . .) with ii+i2 + ■ ■ ■ = ji+ j2 + ■ ■ ■ survive. 



3 Ordinary sum rules and their complements 

Sum rules result from the fermionic representation of the 2d OOP by specific transforma- 
tions of anticommuting variables, keeping the composite nature of the action Sy (21b). 
Let us first rescale by a constant one of the field components, say {■C^^''}; 

z = 0,l,...,iV-l (28) 

Jacobian of the transformation equals to and the fermionic action Sy transforms to 
fiSy- Consequently, 

E 

i,j=0 



Zy = fi-^ fvijVC explfi ^^WiJ'^^ (29a) 

V i,j=0 J 
I- ( 7(JV-1) \ 

= /i^^+Ml?^I?e2.^,exp U Sfc^/ifc^^n (29b) 

\ k,l=0 / 



etc. Zy, a Grassmanian integral, is independent of /i, thus its derivative with respect to 
fi is zero for any value of fi. For the special case /i = 1, the equality din Zy /dfi\^=i = 
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implies 

7{iV-l) 

-N+ Wi,{Ei<i/j) = (30) 

i,j=0 

which, after substituting (p3|), regarding (^5[) and setting = no|V^|, results in the trivial 
neutrality condition 

/ dV [nv{r)-no] = (31) 
Analogously, the equality (9[Zv(Si\E'j)]/(9/i|^=i = yields 

7(7V-1) 

-(iV-l)(S,^,)+ ^h(S.^A^^) = (32) 

fc,/=0 

which is readily shown to be equivalent to the neutrality relation (p!oD . 

Let us now consider another linear transformation of all ^-field components, 

, = o,l,...,iV-l; « = 1,...,7 (33a) 

or all i/j-field components, 

^Pf-^\^^f i = 0,l,...,N -1; a = l,...,7 (33b) 

Jacobian of both transformations equals to A'''^'-^"^-*''^ and the action transforms as Sy — *■ 
Yhj^o^^ ^^'^iWij'^j under (|33a]) and as Sy — > J2]jLo^^ ^-^'^iWij'^j under (|33b| ) . Thus, 

■7(7V-1) 



Zy = A-^^(^-^)/^ /P^^e exp Y A^2,«;,,vl/J 

\ i,j=0 I 

Zy = x-^N{N-i)/2 / p^p^ ^ VZiWi,-^, (34a) 

V jj=0 / 

\ k,i=Q J 
= A-^^(^"i)/2+^yPV^^eS.v|/,exp 5: X'^kWui^A (34b) 
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The equality din Zv/dX\x=i = implies 



1 liN-l) 

--7iV(iV-l)+ Yl ^^..(2.^.) = 

^ id=0 
1 7(A^-1) 

--7iV(iV-l)+ 5: = 

^ ij=0 



(35a) 
(35b) 



On account of (0), this is equivalent to a couple of complex-conjugate equations 
--fN{N-l) + N = Jdhw{z,z) J2 {^i'^j)d^{z'+^z') 



i,j=0 
7(Af--l) 



'yN{N-l) + N = Jd'^zw{z,z) J2 (^i^i)^ 



(36a) 
(36b) 



where we have introduced derivative operators 

1 



2 \dx i (9?/ 



i,j=0 



1 

2 Var 



i 9?/^ 



(37) 



(9"'" = dz,d = dz). They act on complex coordinates according to 

d-^z = 1, d^z = 0; d-z = 0, = 1. 

With \iaw{z, z) = rno0f,(r), 0^ given by (|l^, it is easy to verify the validity of equalities 

1 



d'^\nw(z,z) zuq = / d'^z d \nw(z,z) zuq = — 'jN'^ (3t 
-I jy L -I 2 



d'z 



After some algebra, Eqs.(|3q) then take the form 



1 - 
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iVli-2 



/ d'z 


d'^ In w{z, z) 


z 


Iv 






1 d^z 


d" In w(z, z) 


z 


Iv 







with 5nv{z,z) = nv{z,z) — tiq. Let us denote by dV the positively oriented contour 
enclosing domain V: dV is defined parametrically as follows x = X{(f)),y = F(0);0o < 
< 01. In particular, the ellipse contour admits the parametrization X(0) = acos0, Y = 
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6sin0;O < < 27r. Integrals over V^-domain can be expressed in terms of (9V^-contour 
integrals according to formula 



Jv \ ox oy J Jdv 



(40a) 



where 



dV 



dV 



P{x, y)dx 
Q{x,y)dy 



4'o 



d<j) p[x(</)),y(0)]x'(</.) 

d0Q[X(0),F(0)]r(0) 



(40b) 



Thus, summing and subtracting Eqs.(p9aD and (B9H), one gets respectively 



r^o / d'^r[r- Eb{r)]6nv{r) 
rno / d\[rx Ebir)]Mvir) 



N [2-^^ - d(f) nviX, Y){XY' - X'YlAU) 



d(f) nv{X,Y){XX' + ¥¥') 



(41b) 



where [r x Eb]z = xE^ — yE^. Eq. ( [41a|) was obtained for the case of disk in Ref. 



[see Eq. (4.16)] and represents the generalization of the contact theorem: by removing 
the reference to the boundary and in the radius R oo limit, it results in the contact 
theorem for plane hard wall [0, ||23[. New complementary relation ( [41b|) is informative 
for domain V of generally deformed shape. 



The equality d[Zv{'^i'^j)]/d\\x=i = results in 



---fN{N 



--7iV(iV-l)+j 



7(iV-l) 

k,l=0 
7(7V-1) 

k,l=0 



(42a) 
(42b) 



These relations can be rewritten with the aid of Eqs. (pOD, (|32D and (|35| ) as follows 

7(iV-l) 

{t + l){E,^,) = - J2 {k+l)wki{E,^jEk^if (43a) 

k,l=0 
7(iV-l) 

(j + l)(H,^,) = - {l + l)wM{E.'^jEk'^if (43b) 

k,l=0 
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It is a simple task to pass from ( ^3]) to 
with the obvious generahzation of operators 







= - d^'"2 w 
Jv 






)Z2 


w{ri 


) . 


w{r2) 








= — d^r2 w 
Jv 






)Z2 


win) 


w{r2) 



(44a) 
(44b) 



d^Zj = 5ij, dfzj = 0; zj = 0, zj = 5ij. 



Summing and subtracting (|44a|) and ( [44b|) one finally arrives at 

priQ dVs [r2- Eb{r2)]Sv{ri\r2) = -2nv (n) - n ■ V lUy (n) 
Jv 

- d0 n^[ri; (X, Y)]{XY' - X'Y) (45a) 

J(f)0 

-(n X Vi)^ny(ri) 

+ f^' d0 n^[ri; (X, r)](XX' + YY') {A5h) 



(5no I dVa [rs x Eb{r2)],Sv{ri\r2) 

J 1/ 



Relation ( [45a| ) with ri = was derived for disk in Ref. ||2^ [see Eq. (4.25)]. In the 



R ^ 00 limit of disk, it is related to the dipole sum rule for plane hard wall The 
complementary Eq. ( 45b| ) is original. 



4 New sum rule 

Let us pose the following question: provided that the anticommuting fields under consid- 
eration {^i°'^}Z=i are mapped onto {^i°'\t)}Z=i by a nearest-neighbor transformation 

= a4l\it) + h^tlit), &\t = 0) = (z = 0, 1, . . . , iV - 1) (46) 

with ttN-i = bo = and t being a free parameter, does there exist a choice of the 
coefficients {a,, bi} for which also the composite variables {Sj} ( P2a| ) transform according 
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to a nearest-neighbor scheme 



dt 



a,S,+i(t)+6iS,_i(t), Ei{t = 0) = Ei [t = 0,1, . . . ,^{N - 1)] (47) 



with a^(Ar_i) = 6o = 0? The answer is affirmative []T9[: it can be proven by a direct 
computation that if one chooses 

tti = A{i + 1) bi = B{N - i) (46') 

the consequent {Sj(t)} fulfil differential Eq. (^Tf ) with 

ai = A{i + l) bi = B[-f{N - 1) + 1 - i] (47') 



To determine Jacobian of the mapping, writting formally the solution of (^6]) as 



N-l 

d"H^) = E (48) 

j=0 

there holds 

dci ■ (t) 

= aiCi+ij{t) + bid^ijit), Cij{0) = Sij (49) 
Jacobian equals to det Cij{t)\f~}Q = |c| for each of ^*^°)-components. Its derivative with 



respect to t is given by 



dt 



where Cij{t) is the cofactor of element Cij{t). In combination with Eq. (|49|) , the orthogo- 
nality condition 

N-l 

E (^kjCij = 6ik\c\ (51) 

j=0 

thus leads to d\c\/dt = 0, i.e., |c| = const = 1 for each ^'^"^-component. We conclude that 
Jacobian = 1. For our purpose it is sufficient to consider transformation (46), (|47|) with 
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A = 1 and B = 0; the explicit solution reads 



&\t) = 2,(t)= E rW"^^ (52) 



The insertion of transformation ([5^ ) into the partition function 

7(7V-1) 



C7(JV-1) ^ 

= VijV^expl [^^ + ti^ + l)Ei+^ + 0{e)\wi,^j\ (53) 

I i,j=0 J 

with S^(Ar_i)+i = automatically assumed, and the consequent application of condition 
dlnZv/dt\t=o = lead to 

7(iV-l) 

{t + l)wi,{Ei+i^,) = (54a) 

i,j=0 
7(7V-1) 

5: (j + l)«;,,(S,;Xl/,+i) = (54b) 

where the second formula originates from the t-transformation of {ip^°'^} anticommuting 
fields. The equivalent couple of complex-conjugate equations 

J d^r [d~^ \nw{r)]nv{r) = J d^r 9"'"riy(r) (55a) 
J d^r [d~ lnw{r)]nv{r) = J d'^r d~nv{r) (55b) 

is expressible by using the previously developed formalism as follows 

Tno [ d\ E^{r)nvir) = - f"^' dcp nv{X,Y)Y' (56a) 

Jv J<f>o 

Tno f d\ El{r)nvir) = d0 ny(X, F)X' (56b) 
Jv J4>o 

The t-independence of 



7(JV--1) 



xexp<^ \^k + t{k + l)'~k+i + 0{f)\wki^i\ (57) 



k,l=0 

15 



and similarly of (Sj\E'j(t))t manifests itself at the lowest level as 



y{N-l) 

k,l=0 
7(iV-l) 

k,l=0 



(58a) 
(58b) 



Due to (|54D , the four-correlators (Sj\E'jSfc+i\E'i) in ( ^8a| ) and in ( ^8b|) can be 

substituted by the truncated ones {Ei"^ jEk+i"^ i)'^ and (Hj\l'jSfc\E';+i)^, respectively. Eqs. 
( p8a| ) and ( ^8b| ) are thus expressible as 



w{ri)di 



w{ri) 
w{ri) 



dV2 w{r2)d2 



- / dVa w{r2)d2 
Jv 



^y(n, ^2) 
w(r2) 

nyin, ra) 
ti?(r2) 



(59a) 
(59b) 



These relations can be further simplified to the form 



J^d^r2[dtlnw{r2)]Sv{ri\r2)/q^ = dtnyin) + J^d^r2 d+n^{n, r2) (60a) 



d r2 [82 \nw{r2)\ Sv{ri\r2)/q = Uvin) + / d r2 82 ny{ri, ra) (60b) 



For the elliptic ^-domain of interest it holds lnw{r) = r?T,o(const — ixTyX^ — TrTyy"^) 
Summing and subtracting Eqs. (||0a) and (^b) one then finds respectively 



2nPnoT^ [ d\2 X2Svin\r2) = T^ny(n) + / dVa— n^(ri, rs) (61a) 

Jv OX\ Jv 0X2 

d 



-2n(3noT^ f d'r2y2Sv{n\r2) = ■T—nv{n)+fd'r2—nUn,r2) (61b) 
JV oyi Jv oy2 

To get the diagonal elements of the dielectric susceptibility tensor xv (O) ' o'^^ applies 
jyd?ri xi to (|61a|) and jyd?ri yi to (|61b|) , with the result 



1 



Xv 



Xv 



a + b 



2tt 



2'KTy ah Jo 
1 a + b 1 '•2^ 



2TrTy ah 2Tc'^noq'^ Jo 



d(f) J^d^n xiSv[ri\{X,Y)]cos(f) (62a) 

d0 / dVi y,Sv[n\{X,Y)]sm<p (62b) 
Jv 
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Here, integration per partes was combined with the sum rule ([3T| ) to obtain 
^dVi xi—nvin) = y^dVi— [xinvin)] - no\V\, 



\V\ = Tiab, and analogously for the ^/-component. Eqs. ([6^ ) become simpler in the 
symmetric case of disk of radius R, Xv ~ Xv = Xv, 

XV = --4 ^ / dVi5v^[(n,0i)|(i?,0)]riCOs(0i-0) 



TT R TinQq^ Jo 2tt Jv 
1 1 1 



TT R TinQq^ Jv 
1 1 1 



I dV X Sv[r\{R,Q)] (63) 
dhy Sv[r\{Q,R)] 



TT R TiriQq^ Jv 

where the dependence of S'y[(ri, 0i)|(_R, 0)] on the angle difference |0i — 0| was taken into 
account. Relation (|63| ) can be formally rewritten as 

The final results (|6^)-(|6^) mean an explicit split of Xv i^^o asymptotic l/(27rTy) part 
[see the prediction (^) of macroscopic electrostatics] and the finite-size correction term. 
To show this fact for disk, with regard to the sum rule (p!0D one can write 



j^d^r xSv[r\{R,Q)] = - J^d^r {R - x)Sv[r\{R,0)] (65) 

Moving the reference to the boundary, x' = R — x and y' = y, the integrals on the rhs 
of (^3]) reflect the dipole moment seen from the boundary, which is known to converge 
to a finite value in the thermodynamic limit. We therefore conclude that the correction 



term ^ 1/R. Owing to a slow power- law decay of correlations along a plane wall p5|,|26 
one has to be cautious when identifying the integrals of type ( |65| ) with their asymptotic 
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hard-wall counterparts. Possible vagaries and a check of Eq. (^) are documented in 
Appendix via the exactly solvable 2d OCP at coupling F = 2. 

In conclusion, although the above results (|62D-(|64D were derived strictly for coupling 
constant T = 2*positive integer, it is reasonable to suppose their validity in the whole 
fluid regime. The extension of the treatment to the case of the charged wall and in 
the presence of image forces is straightforward. Potential generalization of the results to 
higher dimensions and to the TCP requires, due to the lack of the fermionic formalism, 
to search for a new method which, on the one hand, reproduces our findings for the 2d 
OCP and, on the other hand, admits a more general applicability like the linear- response 
theory. 
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Appendix 



When F = 2 (7 = 1), the 2d OCP on disk of radius R is exactly solvable 1T§ , [|T7| , |2| . 
The fermionic correlators are given by 

(Si^j) = —6ij, (Si^jHfe^i) = {6ij6ki - Sa6jk) (Al) 

etc., where {wi} (|27| ) are diagonal interaction strengths, 

Wi= / d^r r^'wir) = TT dtfexp(-t) (A2) 
Jv Jo 



written in the units of ttuq = 1. The dielectric susceptibihty tensor is expressible as 

-^Re 1^ dV r'^n{r) + d^ri dV2 ziZ2ny{ri, r2)| (A3) 



Xv 

TTl 



N-1 



It is straightforward to verify the validity of relations 

f dVr^nfr) = (A4a) 
f d'rj d\,-z,z,nUn,r,) = -J^ ^ (A4b) 

Jv Jv ~{ Wi-i 

so that Xv = wn/ (ttNiun-i)- By integration per partes one derives Wi = iwi-i — 
7iN^ exp{—N). Consequently, 

1 N^-^exp{-N) 
XV = (A5) 

71 Wn-1 



On the other hand, 



Since 



\ f d^rx Sv[r\{R,0)] = R n{R) + f dV x n^[r, {R,0)] (A6) 



N~l D2i+1 

R n(R) = w(R) y (A7a) 



d\ xnl[r,{R,0)] = -w{R) y (A7b) 

JV ^ Wi-i 



one arrives at 



(A^ 



Jv Rwn-i 
Inserting ( |A^ ) into (|63D , the exact result ( |A5| ) is recovered (in the units of vrno = l)• 
As ^ oo, the asymptotic form oiwN-i ( |A^ ) can be calculated by the saddle-point 
method in the gaussian approximation: 



vr 



3/2 



wm-1 ~ ^i?A^^"^exp(-A^) 
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(A9) 



Consequently, the quantity ( |A8|) , transcribed according to (^5]), acquires the finite value 



lim / dV (R-x)Sv[r\(R,^ 



7^3/2 



(AlO) 



as was expected. 



One may be tempted to identify formula (|A10D with its obvious plane hard-wall coun- 



terpart. Using the explicit result p5[ for the hard wall localized at x = (plasma appears 
in the half-space x > 0), 



nix] 



no 



IT Jo 



exp 



-{t-xV2f 



n^{xi,X2; \yi - 2/2! 



-nl exp 



-iXi - X2 



1 + 0(t) 
2" 



■dt 



X 



2 roo 6Xp 
TT Jo 



t-ix^+X2)/V2\ -itiy^-y2)V2 

1 + m 



(Alia) 
(Allb) 

dt 



where is the error function 0(t) = (2/y^) /q exp(— 'U^)dM, one obtains 



— dx X 

q'' Jo J-oa 



dy 5(0, x;?/) 



(A12) 



which differs from ( [AlO ) by factor 2. This discrepancy is intuitively associated with the 
slow power-law decay of correlations along the plane wall: an arbitrarily small deformation 
of the boundary towards the circle has a strong impact on this property. The analytical 
(Debye-Hiickel approximation) and numerical (Monte-Carlo simulation) studies of the 



surface charge correlations for finite Coulomb systems were given in Ref. |2^. The 
asymptotic form of these correlations is expected to depend on the shape of the plasma, 
but to be otherwise universal. The exact F = 2 solution for the "soft-wall" version of the 
2d OCP with a quadrupolar field p8[, corresponding to a very large elliptic background, 
supports this finding. 
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